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ABSTRACT 

A recently reported method of predicting heat -transfer coefficients for the extreme 
conditions encountered in the cooling passages of regene rati vely cooled nuclear rocket 
nozzles is used to predict pressures and temperatures at 34 stations in the Phoebus 2A 
nozzle . The predicted pressure and temperatures at the exit of the cooling passages are 
compared with the measured values of EP-IV nuclear tests and found to be in very good 
agreement. Incremental values of coolant temperatures and pressures, coolant passage 
wall temperatures, and heat flux to the coolant are calculated and shoAvn. It is shown 
that a constant Cg of 0. 026 yields values for coolant exit pressure and temperature that 
are in very good agreement with measured values. 
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SUMMARY 

A recently reported method of predicting heat-transfer coefficients for the extreme 
conditions encountered in the cooling passages of regeneratively cooled nuclear rocket 
nozzles is used to predict pressures and temperatures at 34 stations in the Phoebus 2A 
nozzle. The predicted pressime and temperatures at the exit of the cooling passages are 
compared with the measured values of EP-IV nuclear tests and found to be in very good 
agreement. Incremental values of coolant temperatures and pressures, coolant passage 
wall temperatures, and heat flux to the coolant are calculated and shown. A constant 
gas coefficient of 0. 026 yields values for coolant exit pressure and temperature 
that are in very good agreement with measured values. Maximum wall temperature and 

heat flux are in good agreement with those calculated using the conventional varying C . 
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INTRODUCTION 

The extreme conditions encotmtered in regeneratively cooled nuclear rocket nozzles 
produce severe heat -transfer problems in the coolant passages. An effective method of 
predicting heat-transfer coefficients in the cooling passages is essential to the optimiza- 
tion of any nozzle design. Of particular concern is the high heat-flux throat region where 
fluxes of 20 Btu per second per square inch (32. 7 mW/m^ and higher may be reached. 

A munber of experimental investigations have been conducted with single-phase hy- 
drogen flowing turbulently throi^h tubes for a wide range of conditions approximating 
those encountered in the cooling passages of a nuclear rocket nozzle. Each investigation 
resulted in a correlation lor heat-transfer coefficients, which was limited to a particular 
range of conditions. All these investigations are reviewed in reference 1, and a sin^e 



correlation equation for predicting heat-transfer coefficients over a much greater range 
of conditions was reported therein. 

The use of the correlation equation from reference 1 along with the correction fac- 
tors for entrance effects (ref. 2) and curvature (ref. 3) as a method of predicting heat- 
transfer coefficients in the cooling passages of nuclear rocket nozzles was recommended 
in reference 4. An existing digital computer program for calculating heat transfer and 
fluid flow in convectively cooled rocket nozzles (ref . 5) was revised to incorporate the 
recommended heat-transfer correlation equations from reference 4. 

The purpose of this report is to compare the values calculated by the revised com- 
puter program with experimental values obtained from nuclear rocket tests. The 
Phoebus 2A EP-IV nuclear tests were selected because of its wide range of test con- 
ditions. The thermal power of the Phoebus 2A EP-IV tests varied from 490 to 
4080 megawatts. 


PHOEBUS 2A TESTS 

Some overall dimensions of the Phoebus 2 rocket nozzle are shown in figure 1. 
Figure 2 is a photograph of the Phoebus 2A nozzle. The design, fabrication, and nonnu- 
clear testing of the Phoebus 2 nozzles is reported in reference 6. 



Figure 1. - Phoebus 2A nuclear rocket nozzle with some important dimensions. (Ail dimensions in inches (cm).) 
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TABLE I. - OPERATING CONDITIONS FOR PHOEBUS 2A EP-IV NUCLEAR ROCKET TESTS^ 


(a) U. S. Customary Units 


Test 

Flow rate, 
Ibm/sec 

Chamber conditions 

Nozzle inlet conditions 

Thermal power, 
MW 

Hot gas 

Coolant 

Pressure, 

psia 

Temperature, 

°R 

Pressure, 

psia 

Temperature, 

°R 

1 

108.6 

97.9 

124 

1240 

197 

40.8 

490 

2 

142.0 

125. 6 

214 

2120 

362 

43. 0 

1090 

3 

192.2 

176,6 

351 

3080 

650 

47.2 

2190 

4 

244. 1 

219.8 

.435 

2940 

820 

50. 9 

2650 

5 

245. 2 

219.5 

506 

3890 

975 

52. 1 

3630 

6 

262.3 

234. 9 

555 

4060 

1080 

53.9 

4080 


(b) S.L Units 


Test 

Flow rate, 
kg/sec 

Chamber conditions 

Nozzle inlet conditions 

Thermal power, 
MW 

Hot gas 

Coolant 

Pressure, 

MN/m^ 

Temperature, 

K 

Pressure, 

MN/m^ 

Teniperature, 

K 

1 

49.3 

44. 4 

0. 855 

689 

1.36 

22.6 

490 

2 

64.4 

57.0 

1. 48 

1178 

2.50 

23.9 

1090 

3 

87. 2 

80. 1 

2.42 

1711 

4.48 

26. 2 

2190 

4 

110. 7 

99.7 

3.00 

1633 

5. 65 

28.3 

2650 

5 

111.2 

99.6 

3.49 

2161 

6.72 

28 . 9 

3630 

6 

119.0 

106.6 

3 . 83 

2256 

7. 45 

29. 9 

4080 


^Unpublished data received from Donald L. Hanson, Los Alamos Scientific Laboratory. 


Because the Phoebus 2A EP-IV tests covered the widest range of conditions of any 
other Phoebus tests, they were selected for use in this investi^tion. The test conditions 
for the EP-IV tests are shown in table I. 


METHOD OF CALCULATION 


The digital computer program reported in reference 5 was revised to incorporate 
the recommended heat-transfer equation from reference 4. 









































Heat-Transfer Calculations 


Heat-transfer equations for both the hot-gas side and the coolant side are required to 
obtain a heat balance through the coolant wall. It is important that the best available 
prediction equations be used on both sides of the coolant passage wall. The correlation 
equation for the hot-gas side has not been improved upon but the coolant side correlation 
has undergone considerable improvement (refs. 1 and 4). These changes will be dis- 
cussed in the section on the coolant side, 

Hot- gas side. - The Nusslet equation 


NUf = CgRe^- ®Pr^- 3 (1) 

is used in the computer program. Common practice (ref. 5) is to use a that is a 

function of nozzle area ratio. In this investigation both a constant G of 0.026 and a 

& 

TABLE n. - AREA RATIO AND GAS COEFFICIENT 


FOR EACH CALCULATION STATION^ 


Region 

Station 

Area 

Gas 

Region 

Station 

Area 

Gas 



ratio 

coefficient, 



ratio 

coefficient, 




C 










g 





g 

Divergent 

1 

7.00 

0.032 

Convergent 

21 

2. 34 

0.027 


2 

6. 44 

.033 


22 

3. 

17 

. 028 


3 

5. 90 

. 033 


23 

4. 

13 

.028 


4 

5.39 

.033 


24 

5. 

21 

.029 


5 

4.90 

.034 


25 

6.42 

.032 


6 

4.43 




26 

7. 

76 

.035 


7 

3.98 




27 

9. 22 

. 038 


8 

3 . 56 









9 

3. 17 



Knuckle 

28 

10.81 

0. 042 


10 

2.79 




29 

11. 85 

.045 


11 

2.44 




30 

12. 22 

. 048 


12 

2. 12 

.032 

Chamber 

31 

12. 22 

0. 056 


13 

1.81 

. 030 


32 



. 064 


14 

1.53 

. 028 


33 



.075 


15 

1.28 

. 024 


34 



.080 

Throat 

16 

1.04 

0. 018 

1 






17 

1.00 

.018 







18 

1.07 

.028 







19 

1.27 

.029 






i 

20 

1 

1.63 

.029 







^Unpublished data received from James O. Sane, Aerojet-General Corp. 
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variable Cg are used and the results compared. The variable Cg and the area ratio 
for the 34 stations in the nozzle is shown in table II. The equation for the gas-side heat- 
transfer coefficient used in the computer program is 


C„K f 

h = ^ 

S n 

^H,g 


^w\%Pf 

[AfJ 


o- ^ ^ 


The heat-transfer surface area used in the computer program is 


( 2 ) 


Aht AiTTRj.^ av^c, av 


(3) 


where is a heat-transfer area cprrection, which compensates for the variation of 

heat flux around the perimeter of the coolant tube. The generally accepted of 

value 0. 8 (ref. 5) is used in this investigation. ’ 

Although the radiative heat flux at each station may be inserted as input data to the 
computer program, it was neglected in this investigation. There is no provision in the 
computer program for the small amount of film cooling on the hot-gas side of the 
chamber wall. The maximum effect can be estimated by assuming no heat addition at 
the last 4 stations (about 12. 5 in. or 32 cm). 

Hot- gas static temperature and pressure calculations for hydrogen at equilibrium 
conditions for assigned chamber pressure, chamber temperature, and nozzle area ratios 
were made using the computer program of reference 7. 

Coolant side. - A Nusselt type correlation equation with the physical properties and 
density evaluated at the film temperature is part of the computer program in reference 5. 
This correlation equation was replaced by the one reported in reference 1 and used as 
recommended by reference 4. The correlation equation from reference 1 is 



0.023 Re»;fpr«;j% 


(4) 


where 



57-1.59/(Z/Djj^p] 


G 


(5) 



The correlation for entrance effects from references 2 and 4 is 





where F = 11 for an orifice entrance. The correction for curvature effects from ref- 
erences 3 and 4 is 


C 


5 ■ 



(7) 


where Cg = 0. 05 for the concave or swept surface and -0. 05 for the convex or unswept 
surface. Equations (4) to (7) are combined to give the single correlation equation 



= 0. 023 Re°- ^Pr°’ fc^C.C. 

b,l ^ ^ ^ 


( 8 ) 


which is used in the computer program with and Cg being used only where appro- 
priate. Figure 3 shows the form of equation (8) used in the various coolant flow regions. 



NU[,_j= 0.023 Re§;?Pr§;fC 3 C 4 C 5 Region 



1 €5 = 1 

2 C4=1,C6=0.05 

3 C4» l, Cg ' l 

4 C4- 1, Cg- -0.05 

5 C4- 1, 05 -1 






H.i 


Figure 3. - Five regions of the coolant passage and 


constants applicable to each region. 
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Fluid Flow Calculations 

Both friction and momentum pressure drops were calculated in the computer pro- 
gram. The eqiiations used to calculate these pressure drops in reference 5 were not 
changed and are repeated here for the reader’s convenience. 

Friction presswe drop. - The friction pressure drop was calculated usii^ the equa- 
tions already in the computer program reported in reference 5. The friction pressure 
drop was calculated by the equation 


Ap,^ = 


^s, av^H, av^ 


(9) 


Two equations for calculating friction coefficients are available in the program: one 
for smooth tube conditions, the other for rough tube conditions. Since the Phoebus 
nozzle cooling passages have a relative roughness e of 60 microinches (1. 52 pm), the 
rough tube equation 


-L=-4,0Iog 2 + 1-255 (JJ) 


is used. 

The friction coefficient is increased by the curvature of the throat and knuckle re- 
gions. The Ito correction factor (ref. 3) for curvature 


C 


7" 



( 11 ) 


is applied as follows: 


with f being used in equation (9) . 

Momentum pressure drop . - The momentum pressure drop eqiiation reported in 
reference 5 is 


8 




( 12 ) 


It is used in this investigation to calculate the momentum pressure drop in the coolant 
passages of Phoebus 2A nozzle. 


DISCUSSION OF RESULTS 

Coolant temperatures and pressures, wall temperatures, and heat fluxes at 34 sta- 
tions along the coolant passages were calculated using the revised computer program 


Film coolant 



34 32 28 24 20 16 12 8 4 0 


Station 

Figu re 4. - Variation of coolant total tempera- 
ture and static pressure with axial location. 
Thermal power 4080 megawatts; chamber pres- 
sure to 555 psia (3. 83 MN/m^); chamber tem- 
perature, 4060° R (2256 K). 
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from reference 5, Total temperature and static pressure of the coolant at each axial 
station for test EP-IV-6 are shown in figure 4. The only measured values for compari- 
son with the calculated temperatures and pressures are taken at the exit of the coolant 
passages (the inlet of the reactor core reflector). The measured values for test EP-IV-6 
are shown in figure 4, The solid line represents the calculated values with Cg varying 
as shown in table II. The dashed line represents the calculated values using a constant 
Cg of 0. 026, The agreement between both calculated values and the measured values is 
very good. The measured and calculated exit pressures and temperatures and percent 
deviation for several other EP-IV tests are shown in table III. For the case of no heat 
addition in the last 4 stations (corresponding to film cooling of the chamber wall), both 
the exit temperature and pressure fell within the range of accuracy of the measured 
values and are not shown in table m. 

In figure 5 the local hot-gas side wall temperature and the local coolant side heat 
flux is shown as a function of axial position. The solid line represents the values calcu- 
lated using the variable Cg from table n and the dashed line represents the values re- 
sulting from the use of a constant Cg of 0. 026. The greatest differences resulting from 
the calculations using the different Cg’s are in the chamber. The use of temperature 
rise and pressure drop measurements indicates little of what happens locally in a rocket 
nozzle. 

The maximum wall temperature is 1936'^ R (1076 K) for a constant Cg of 0. 026 
compared with 1838° R (1021 K) for the variable Cg from table II with both appearing at 
approximately the same location. The maximum heat flux to the coolant is 17. 4 Btu per 

second per square inch (28. 4 MW/m ) for C = 0. 026 and 16. 9 Btu per second per 

o & 

square inch (27. 6 MW/m ) for Cg from table II. The maximum heat flux and maximum 
wall temperature do not appear at the same axial location, 

Calculations using either a variable C or a constant Cg of 0. 026 both give coolant 
exit pressures and temperatures that are in good agreement with measured values. The 
maximum heat flux and wall temperature for the constant Cg varies little from those 

computed with a variable C . Since the variation of C_ with area ratio varies a great 

& & 

deal with both investigator and nozzle configuration (ref. 5), the capability of using a 
constant Cg for reliable nozzle calculations is of great value. 

The rather abrupt changes in wall temperature is due to the application of the curva- 
ture correction which increased the heat-transfer coefficient on the concave or swept 
surface (decreasing the wall temperature) and decreased the heat-transfer coefficient on 
the convex or unswept surface (increasing the wall temperatxire) . In the actiml nozzle 
wall these changes would probably not be so abrupt because of axial heat conduction in 
the coolant passage wall. Axial heat conduction was not accounted for in the computer 
program. 
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TABLE III. - CALCULATED AND MEASURED TEMPERATURES AND PRESSURES 



Eh 




ci • 
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Unpublished data received from Donald L. Hanson, Los Alamos Scientific Laboratory. 










Film coolant 




Figure 5. - Variation of hot -gas side wall temper- 
ature and coolant side heat flux with axial loca- 
tion. Thermal power, 4080 megawatts chamber 
pressure, 555 psia 1^83 MN/m^); chamber 
temperature, 4060° R (2256 K). 


CONCLUDING REMARKS 

Recently reported equations, which accurately predict coolant heat-transfer coeffi- 
cients (refs. 1 and 4), have been inserted into an existing computer program (ref. 5) to 
calculate coolant exit temperature and pressure. These calculated values are in very 
good agreement with measured values. One set of calculations used recommended G^’s 
that varied from 0. 018 to 0. 080, another set of calculations used a constant Gg of 
0. 026. Both sets of calculations yielded exit temperatures that are in good agreement 
Avith measured temperatures. The constant calculations predicted exit pressures 
that were in a little better agreement with measured pressures than the variable Cg 
calculations. These calculations show that with accurate predictions on the coolant side 
the less well developed heat-transfer correlations on the hot-gas side are of secondary 
importance in predicting total temperature rise and static pressure drop. The hot-gas 
side heat-transfer does change the wall temperature distribution, which indicates that 
measured wall temperatures are needed to verify any hot-gas side equations. 
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The close agreement between predicted and actual performance gives confidence 
that now these prediction equations can be used to design nuclear rocket nozzles that are 
more reliable and possibly lifter in weight. 

The small effect of whether a variable or a constant Cg was used in the calculation 
of coolant exit pressure and temperature and maximum heat flux and gas -side wall tem- 
perature indicates that reliable calculations can be made using a constant Cg of 0. 026. 
Since C has always varied with both investigator and nozzle configuration, the capa- 

o 

bility of using a constant C greatly simplifies nozzle calculations. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, ]\fey 27, 1969, 

122-28-02-33-22. 
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APPENDIX - SYMBOLS 

area 

ntimerieal coefficient 

convective heat-transfer coefficient correction for fluid properties variation 
convective heat-transfer coefficient correction for entrance effects 
convective heat-transfer coefficient correction for curvature effects 
exponent of Cg, 0. 05 for concave surface, -0. 05 for convex surface 
friction coefficient correction for tube curvature 
specific heat at constant pressure 
diameter 

relative roughness of surface 
entrance effect coefficient 
friction factor for straight tubes 
friction factor for curved tubes 
mass flow per unit cross-sectional area 
gravitational conversion factor 
convective heat-transfer coefficient 
thermal conductivity 

linear distance along coolant passage wall 

linear distance along coolant passage wall between stations 

Nusselt number 

Prandtl number 

pressure 

friction static-pressure drop 
momentum static-pressure drop 
local heat flow rate 
radius 

radius of curvature 
Reynolds number 
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T temperature 

w mass weight of flow 

€ heat 'transfer surface-area correction factor 

c 

fi dynamic viscosity 

p densiiy 

Subscripts: 
av average 

b bulk 

c calculated 

f film 

fl flow 

g gas 

H hydraulic 

ht heat transfer 

i inside 

I liquid or coolant 

m measured 

o outside 

s static 

t tube 

w wall 
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